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The Glauber approximation for the treatment of heavy-ion scattering, 
has already been shown to give reliable predictions for the reaction cross sec- 



X 



tion in the particular case of intermediate energy charge-exchange processes. 



In the present work, we couple a Glauber-type model to microscopic Ran- 
dom Phase Approximation calculations of the charge-exchange excitations of 
208 Pb. The aim is to solve the longstanding question whether the very elu- 
sive charge-exchange isovector monopole has been really identified in the past 

experiments, or other multipoles were prevalent in the observed spectra. 
PACS numbers: 21.60. Jz, 24.30.Cz, 25.55.Kr 
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I. INTRODUCTION 



While some systematics is available for isoscalar giant resonances in atomic nuclei, many 
isovector charge-exchange modes still remain elusive. A number of intermediate energy 
(i.e., at bombarding energies Ei a b larger that 100 MeV/A) reactions favour the excitation 
of spin-flip modes due to the dominance of the a • a part of the effective nucleon-nucleon 
interaction. For selective excitation of the non spin-flip modes, either pions or low-energy 
reactions or special projectile-ejectile combinations are needed. Pioneering experiments of 
the first two types (cf., respectively, Refs. [1] and [2]) showed resonance-like structures, 
but the multipolarity assignments were difficult or impossible. This contrasts with the 
abundance of theoretical predictions for the charge-exchange giant resonances obtained in 
particular within the framework of self-consistent Random Phase Approximation (RPA), 
starting from the work of N. Auerbach and A. Klein [3]. The lack of knowledge contrasts 
as well with the importance of charge-exchange isovector resonances for many purposes. 
The isovector giant monopole resonance (IVGMR) is related to the isospin impurity of the 
nuclear ground state [4] , and to the isovector nuclear incompressibility [5] - a nuclear matter 
parameter which is essentially unknown. More generally, the study of the excitations which 
involve the isospin degree of freedom, should give the possibility to test the nuclear models 
and eventually improve their predictive power for extreme conditions (like nuclei with large 
isospin values, or nuclear matter close to isospin instability). 

The most recent experiments aimed to the identification of the IVGMR and/or other 
multipole resonances, have employed the ( 13 C, 13 N) reaction. This reaction is believed to 
have good selectivity for the non spin-flip type of excitation in the target, because of the 
dominance of Fermi over Gamow- Teller transition in 13 C. In the first experiment performed 
in GANIL [6] at 50 MeV/A, the candidate for the monopole excitation in 208 T1 (this is the 
AT Z =+1 excitation of the 208 Pb target) showed up at 21 MeV 1 , in large disagreement with 



The excitation energies quoted in this work, are referred to the ground state of 208 Pb, which 
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the result of 11 MeV found in the pion experiment of Ref. [1]. One has to remember the 
severe problems of this pion experiment (low statistics, large background, full width of the 
monopole peak comparable to its energy), but also in the ( 13 C, 13 N) experiment there was the 
puzzling feature of the angular distribution which was inconsistent with AL=0. In a second 
experiment [7] at 60 MeV/A, the neutron decay of the monopole candidate at 21 MeV was 
measured, with the hope that an isotropic angular distribution for the emitted neutrons 
could support the multipolarity assignment. However, the results from this measurement 
are still ambigous, as we discuss in more detail in Sec. IV. 

Concerning the comparison with theory, the calculations of Refs. [3,8] reveal a marked 
discrepancy of self-consistent RPA calculations with the results of the ( 13 C, 13 N) experiment. 
According to Ref. [8], the monopole peak is sensitive to the effective force employed (a 
Skyrme parametrization), and it is found either at 13 or 16 MeV, therefore much lower than 
the energy region in which the experimental candidate for IVGMR is found. On the other 
hand, in the same paper it was pointed out that in that energy region the RPA calculations 
predict a sizeable amount of octupole strength. A calculation of the relative values of AL=0 
and AL=3 cross sections is therefore highly desirable. This further step of calculating the 
cross sections (instead of simply limiting the analysis to the strength distributions) finds its 
justification in the fact that the former do not obviously scale as the corresponding strengths. 
Simple reaction frameworks that can however directly exploit the structure information 
provided by the transition densities are available in the literature, as for example the one 
offered by the Glauber model. As shown in Ref. [9], the extension of the Glauber model 
to the the description of heavy-ion charge-exchange reactions at intermediate energy, like 
in the case at hand, can be rather successful in terms of comparison with experimental 
data. This motivates the present calculation of the cross sections for electric multipole 
charge-exchange excitations in the Glauber approach, at 60 MeV/A, having as ingredients 
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microscopic transition densities of the RPA self-consistent calculation. 

We believe the present work may shed light on the experimental problem of the IVGMR 
identification. We should mention that the efforts to improve our rather poor knowledge of 
the various charge-exchange modes still continue in different laboratories [10]. 

The outline of the paper is the following. The formalism for the extension of the Glauber 
model to heavy-ion charge exchange reactions is briefly recalled in Sec. II. The RPA calcu- 
lations for the different multipolarities are presented in Sec. Ill, where they are discussed in 
particular in connection to the coupling to the continuum and to the two particle-two hole 
states. The predictions of the model for the charge-exchange cross sections are presented 
and discussed in Sec. IV, with special emphasis on the energy range selectively populated 
by each multipolarity. The main conclusions of the paper, as well as the implications for 
future experimental and theoretical work, are finally given. 

II. CHARGE EXCHANGE CROSS SECTIONS 

High energy scattering processes are dominated by nucleon-nucleon collisions and can 
therefore be well reproduced by the Glauber model [11]. In this approximation, nucleus- 
nucleus reactions are fully microscopically described in terms of the nucleon-nucleon scat- 
tering amplitudes and nuclear (transition) densities. The model was developed in detail 
for charge-exchange reactions in a recent paper [9] where the formalism is presented and 
is applied to jo-shell nuclei. For the sake of completeness we recall here some of the main 
features and formulae relevant for the present work. 

We assume that at relatively high energies these reactions are dominated by one step 
processes [12]. Within this approximation, the scattering amplitude of one-step charge- 
exchange reactions (Aa — > Bb), is expressed as, 

f Aa ^ Bb (A) =ikJ2[bdb ^lT Bb (b) J M (Ab) e^+^ e -iM* A> (21) 

LM J 

where A is the transferred momentum, k is the nucleus-nucleus relative momentum, b is the 
impact parameter, X(b) +x(6) represents the phase shifts including the nuclear and Coulomb 
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contributions and Jm(A6) is the Bessel function of order M. The matrix element ntlT Bb (b) 
takes into account the details of the process and in Ref. [9] a general expression for any type 
of charge-exchange reactions was given. Since in the present case we will be only considering 
non-spin-flip transitions (AS=0), we cast the expression of plm (Eq. (2) of Ref. [9]) in the 
form, 

ni a M Bb { h ) = ~7~ — ( tan a t b -n b | 1 n a - n b )( t A n A t B - n B \ ln A -n B ) 



ik 



NN 

E ( 3a m a j b -m b | J p m a - m b )( j A m A j B - m B \ J t m A - m B ) 
JpJt 

[Bj v B Jt ] L M J q dq p J a l{q) p J Ud) JnAq) M^), (2.2) 

where kNN is the nucleon-nucleon relative momentum. The above expression includes the 
scattering amplitude f$ N , related to the isovector central part of the effective nucleon- 
nucleon potential and the Fourier transforms p(q), of the projectile and target transi- 
tion densities p'ib( r ) an d P AB ( r ), respectively. The quantities B arise from the projec- 
tion of the spherical harmonics on the plane perpendicular to the trajectory (Y^ m (|,<^) = 
i t+m Bt m e imip 8t +mteven ). The squared brackets are used to indicate the coupling of the asso- 
ciated angular momenta. 

From the Eqs. (2.1) and (2.2) we can straightforwardly obtain the differential cross section 
to each final state, 

%(Aa^B b ,0) = 
1 ( k" 2 



, ( t a n a t b - n b | 1 n a - n b ) 2 ( t A n A t B - n B \ 1 n A - n B f 
J2 E [B.J p Bj t ] L M [ b db e x ^ + ^ J M (Ab) e- m ^ [ q dq p J » p Jt f nn (qb)J M (qb) 

J p Jt LM J 

and, integrating over the scattering angle, the corresponding total cross section 



(2.3) 



VAor+m = J bdb \e^ + ^ pitT Bb (b)\ 2 . (2.4) 



The nucleon-nucleon scattering amplitudes have been taken from the tabulation given 
in Ref. [13] and interpolated for the energy of interest, i.e., 60 MeV/A. The transition 



densities corresponding to projectile (ejectile) have been obtained as in Ref. [9] using the 
wave functions that are solution of a shell model calculation in the p-shell [14]. Due to the 
non spin-flip character of the processes involved, only transitions with J p = are allowed 
for the ( 13 C, 13 N) reaction regardless of the states populated in 208 T1. The calculation of the 
other crucial ingredients, i.e., the transition densities to each excited state in the target, is 
discussed in the next Section. 

III. MICROSCOPIC TRANSITION DENSITIES FOR TARGET EXCITATION 

As mentioned already, the transition densities corresponding to 208 Pb excitations are 
derived from charge-exchange RPA using Skyrme effective forces. The method is well known, 
and the results for the IVGMR in 208 T1 have been published in Ref. [8]. We recall however 
the main features of this type of calculations and we report some detail of the present one. 

We choose the Skyrme parametrization SIII [15]. In Ref. [8], the sensitivity to the choice 
of the Skyrme parametrization is discussed in the IVGMR case. We come back to this point 
in the next Section. In our method, we first solve the Hartree-Fock (HF) equations in real 
space and we obtain the mean field. The unoccupied states, including those at positive 
energies, are obtained by diagonalizing the mean field on a harmonic oscillator basis (in this 
case huo = 6.2 MeV), that is, the continuum is discretized. Using these single-particle states, 
a basis of proton particle-neutron hole plus neutron particle-proton hole configurations is 
built, and the RPA matrix equations are solved on this basis (which is large enough to ensure 
that the appropriate sum rules [3] are satisfied). Our charge-exchange RPA is self-consistent, 
in the sense that the residual interaction between p-h configurations is derived properly from 
the Skyrme force. Only the spin-orbit residual two-body force is dropped. The number of 
configurations used in the solution for the different multipolarities (that is, AJ=AL=0,1,2 
and 3) is provided in Table I. It corresponds of course to the number of RPA states: all 
these states have been included in the calculation of the cross sections, that is, no selection 
based on their energy and/or strength has been applied. We present the results associated 
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with the excited states of 208 T1. In this case, the states have good isospin T + 1 (T is the 
g.s. isospin of 208 Pb) and it is possible to use them within the framework of our Glauber 
model whose formulation includes the isospin quantum number. 

Solving the RPA gives of course the energies E n of the excited states \n), as well as 
their wave functions in terms of the and amplitudes. The corresponding radial 
transition densities (mentioned above after Eq. (2.2)) are given by 

^ t (0=E(^S ) -^ B) )<Plln||/»>i2p(r)i2 h (r) J (3.1) 

ph 

where ph labels the particle-hole components of the basis and R(r) is the radial part of the 
single-particle wave function. These transition densities are related to the strength of each 
state S n associated with the multipole operators O = Y^irfY L (fi)t^ by 

S n = \ j drr L+2 p J n *{r)\ 2 . (3.2) 

The quantity p% (r) must be multiplied by a factor yj 1 2 ^j^ in order to be consistent with 
Eqs. (2.2) and (2.3). 

In Ref. [8] more sophisticated calculations were performed. Using the discrete RPA as a 
starting point, it is possible to take into account the coupling with the continuum as well as 
with configurations of two particle-two hole type using the model first introduced in Ref. [16]. 
Including the continuum coupling is equivalent to solve the well known continuum-RPA. 
Both couplings have as a consequence a shift of the RPA peak(s), and also a redistribution 
of the strength whose main effect is to provide the escape and spreading widths to the 
RPA peak(s). In the case of the IVGMR in 208 T1, as discussed in great detail in Sec. Ill 
of Ref. [8], the discrete RPA peak is shifted downwards by about 2.8 MeV, from 19 MeV 
to 16.2 MeV. This effect is not included in the present calculation, therefore the peaks 
corresponding to the + cross section in 208 T1 should be moved towards lower energies. This 
reinforces the final message of the present paper, namely that the octupole cross section 
dominates in the region around 21 MeV. In fact, we have checked that the octupole peak 
is much less shifted by the mentioned couplings. To substantiate this point, in Fig. 1 we 
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show the result of a full calculation of the octupole strength, performed exactly as in [8] 
including escape and spreading, compared with the averaging of the discrete RPA results 
with Lorentzian distributions having 2 MeV width. The full calculation gives essentially the 
same distribution as the discrete RPA + averaging. 

IV. RESULTS FOR THE CROSS SECTIONS AND CONCLUSIONS 

The calculations described in previous Sec. Ill provide the radial transition densities to all 
excited + , 1~, 2 + , and 3~ states in 208 T1. As indicated in Sect. II these transition densities 
provide the basic ingredients for the calculation of the charge-exchange cross sections. As an 
application, we discuss here the results of our model for charge-exchange reaction induced 
by 13 C projectiles at 60 MeV/A on the 208 Pb target. 

A first question to be posed, is whether the profile of the cross sections as a function 
of the excitation energy follows that of the multipole strengths. In Fig. 2(a) we show some 
examples that illustrate this point in the case of dipole states. We find that in some cases 
states with similar dipole strength have rather different total cross sections (as for example 
the states at excitation energies E n = 25.5 and 27.3 MeV) and conversely, states with very 
different strength have quite similar cross sections (as states at E n = 27.3 and 29.1 MeV). 
This fact can be understood by inspecting the radial behaviour of the transition densities 
at the nuclear surface [Fig. 2(b)]. While the dipole transition matrix elements involve the 
transition densities in the whole radial space, charge-exchange processes induced by heavy 
ions have a rather grazing character and the reaction amplitudes get their contributions 
mainly from the surface region. We see therefore that the states corresponding to the 
largest cross sections are those which have largest values for their transition densities at the 
surface, even if the associated strengths are not necessarily large. 

On the basis of the above conclusion it is worth studying, as a function of energy, the 
population of the different states produced in a charge-exchange process, rather than to limit 
oneself to the strength distributions. We have calculated the angle integrated cross sections 
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of Eq. (2.4) for all the above indicated multipoles. In order to take into account the widths 
of the states we associate to each discrete RPA state a Lorentzian averaging distribution as 
described in Sec. Ill, accounting for both escape and spreading mechanisms. The widths of 
the Lorentzian distributions were chosen to be 2 MeV. The results are displayed in Fig. 3. 
From the figure one can see that the main contributions correspond to the quadrupole and 
octupole transitions in the energy range below 25 MeV. 

We finally discuss our predictions in more detail in connection with the experiment of 
Ref . [7] . The aim of the experiment was to infer the multipolarity of the observed structure 
from the angular distribution of the neutron decay. From the excitation energy spectra at 
small angles - which in principle favour low multipole transitions, in particular AL=0 - a 
large peak around 21 MeV is observed. The authors themselves point out that the statistics 
for the neutron detection was too low to obtain a clear signature of the multipolarity of this 
resonance from the neutron decay. Nevertheless, they attributed this peak to the presence of 
the isovector giant monopole resonance. From the present analysis we find that the octupole 
excitation should be dominant in the 21 MeV energy region, as shown in Fig. 3 for the total 
cross sections. This is in agreement with the remark made in [8] on the simple basis of 
the strength distribution: nevertheless, we have seen that cross section distributions and 
strength distributions can be different. 

For small angles the relative importance of monopole transitions increases. Therefore, a 
valid question is how the different multipolarities would compare when only forward angles 
are considered. In Ref. [7] two different spectra are shown: one for the cross section inte- 
grated between and 1 degree and another integrated between 1 and 2.5 degrees. Hence, 
we have calculated the excitation functions for the monopole and octupole transitions by 
integrating the cross sections of Eq. (2.3) in the same angular range. The results are shown 
in Fig. 4. One can see that the relative importance of the monopole with respect to the 
octupole excitations is larger at small angles, but the octupole overcomes the monopole 
transition in both angular ranges considered in [7]. 

In conclusion, the present paper shows the feasibility and usefulness of fully microscopic 
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calculations which couple structure and reaction models. The Glauber model, which was 
already shown to be quite reliable for the study of intermediate energy charge-exchange 
reactions, has been applied to the case of 60 MeV/A ( 13 C, 13 N) scattering on 208 Pb. This 
kind of experiment was peformed with the aim of extracting the properties of the IVGMR. 
We have shown that the octupole cross section dominates over the monopole in the energy 
range around 21 MeV, where the experimental peak was identified, if one considers the angle- 
integrated total cross section. The ratio between the monopole and octupole cross sections 
grows if they are both integrated over small angular ranges; nevertheless, the monopole 
is not dominating in the angular ranges considered in the experiment of Ref. [7]. These 
remarks may be of some importance when the quest for the IVGMR, as well as for the other 
still poorly known charge-exchange multipoles, continues. 
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FIG. 1. Distribution of charge-exchange octupole strength in 208 T1. A complete calculation 
including escape and spreading (full line) is compared with the averaging of the RPA results with 
Lorentzian functions having 2 MeV width (dashed line). 
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FIG. 2. (a) comparison of the relative contributions, denned in terms of percentages of the 
total, to the dipole strength (light shadowed bars) and to the cross section (dark shadowed bars) 
in the case of three different 1~ states in 208 T1, whose energies are reported in the figure, (b) 
Behavior of the corresponding transition densities at the nuclear surface. 



13 



6x10- 2 




10 20 30 40 50 



E (MeV) 



FIG. 3. Total cross sections as a function of the energy in the residual nucleus 208 T1 for the 
different multipolarities. 
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FIG. 4. Cross sections as a function of the excitation energy for monopole and octupole tran- 
sitions in 208 T1 integrated in the angular range: (a) to 1 degree, and (b) 1 to 2.5 degrees. 
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TABLES 

TABLE I. Number of particle-hole (ph) states included in the RPA calculation of the different 
multipole excitations. 



J n neutron p-proton h proton p-neutron-h 
0+ 64 88 
1- 172 240 
2+ 244 348 
3~ 288 420 
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